REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway, 

Suite  1 204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

29  November  2016  Briefing  Charts 

3.  DATES  COVERED  (From  -  To) 

02  November  2016  -  30  November  2016 

4.  TITLE  AND  SUBTITLE 

A  Hybrid  Model  for  Multiscale  Laser  Plasma  Simulations  with  Detailed 
Collisional  Physics 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

David  Bilyeu,  Carl  Lederman,  Richard  Abrantes,  Hai  Le 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

Q02Z 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/RQRS 

1  Ara  Drive 

Edwards  AFB,  CA  93524-7013 

8.  PERFORMING  ORGANIZATION  REPORT  NO. 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/RQR 

5  Pollux  Drive 

Edwards  AFB,  CA  93524-7048 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 

NUMBER(S) 

AFRL-RQ-ED- V  G-2016-328 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  Unlimited.  The  U.S.  Government  is  joint  author  of  the  work  and  has  the  right  to  use, 
modify,  reproduce,  release,  perform,  display,  or  disclose  the  work. 


13.  SUPPLEMENTARY  NOTES 

For  presentation  at  AFOSR  Plasma  and  Electroenergetics  Review  Meeting;  Arlington,  VA  (29  November  2016) 
PA  Case  Number:  #16531;  Clearance  Date:  07  November  2016 
Prepared  in  collaboration  with  ERC 

14.  ABSTRACT 

Viewgraph/Briefing  Charts 


15.  SUBJECT  TERMS 

N/A 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE 

PERSON 

D.  Bilyeu 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

SAR 

23 

19b.  TELEPHONE  NO  (include  area 
code) 

N/A 

Standard  Form  298  (Rev.  8- 
98) 

Prescribed  by  ANSI  Std.  239.18 


A  Hybrid  Model  for  Multiscale  Laser  Plasma  Simulations  with 

Detailed  Collisional  Physics 


AFOSR  Plasma  and  Electroenergetics  Review  Meeting 

29-30  November  2016 


Richard  Abrantes1’2,  Hai  Le1’2, 

Carl  Lederman2, 

Collaborators:  PI"  David  Bilyeu3 

’UCLA 
2ERC  Inc. 
3AFRL/RQRS 


•  R.  Caflisch  (UCLA) 

*  A.  Karagozian  (UCLA) 


Distribution  A  -  Approved  for  public  release;  distribution  is  unlimited;  PA:  #16531 


Outline 


*£, 


*  Goals 

*  Review  of  Past  Work 

*  Argon  Collisional-Radiative  Complexity  Reduction  Validation 

*  Non-Maxwellian  CR 

*  Phase-accurate  Multiscale  Particle  Push 
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Goals 


Utilize  hybridization  techniques  to  produce  accurate  and  efficient  plasma 
simulations  that  spans  many  orders  of  magnitude  in  both  space  and  time. 

Capture  complex  physics:  excitation/ionization,  transport,  radiation,  etc. 

Consistent  collision  operator  across  different  levels  of  fidelity. 


Current  Focus: 

•  Generalization  of  collisional- 
radiative  kinetics  with  level 
grouping 

•  General  Hybridization  techniques 

•  Focus  on  each  solver  before 
hybridization 

•  Special  attention  to  low  density 
low  energy  conditions 


Laser  Plasma  Interaction 


FRC 
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Summary  of  Past  Work 


Maxwellian  Inelastic  Collisions 


•  Detailed  CR  model  for  multiple  ionization  stages 

•  Validation  against  experimental  data 

•  Nonequilibrium  radiation  transport:  coupling  with  a 
collisional-radiative  model 

•  Inelastic  collisions  in  a  MF  plasma:  enhanced 
thermochemical  kinetics. 


Multiscale  Hybridization 


•  A  time-parallel/multiscale  method  with  energy 
preservation 


Analytical  BGK 


•  Conservative  even  for  large  mass  ratios 

•  Conservation  are  independent  of  collisional 
frequencies 
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Coll  sional  Radiative  (CR)  Overview 


Updates 


•  Expanded  complexity  reduction  to  include 
multiple  ionization  levels 

•  Adaptive  integration  technique  ->  fixes  rate 
calculation  for  higher  electron  temperature 

•  Investigate  grouping  sensitivity 

•  Linked  with  LANL  database  for  Argon  cross 
sections  and  atomic  level  information 

•  Algorithms  not  hard  coded  for  Argon. 


Levels  of  Complexity 


•  Full  rates  (LANL) 

•  Cutoff  nearly  ionized  levels 

•  Grouping  Strategies 

•  Uniform 

•  Boltzmann 

•  QSS  (Boltzmann  and  Planck  equilibrium) 

•  Group  Selection 

•  Electron  configuration  (no  splitting 
information) 

•  Highly  excited  states 

•  Analysis  of  full  run 

•  Numerical  optimization 
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CR  Governing  Equations 
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CR  Governing  Equations  cont. 


•  Current  Assumptions 

•  Electron  dominated  collisions 

•  6  function  for  ion  distribution 

•  Maxwellian  electrons 

•  Current  Model  Includes 

•  Multiple  ionization  levels 

•  Excitation/de-excitation 

•  lonization/recombination 

•  Multi-photon  ionization  and 
inverse  Bermsstrahlung 

•  Radiation  losses  via,  Bound- 
Bound  and  Bound-Free 


*previous  included 
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CR  Grouping  Techniques 


Uniform 


•  Solely  conserves  number 
density 

•  Weights  each  level 
according  to  level 
degeneracy  within  group 


Conserved  Variable: 

Nn  =  YjNi 

ien 

Effective  rate: 
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dN 


n 


dt  Ne 


^  a(m\n)Nn  +  ^  P(m\rC)Nn 


m>n 


m<n 


+ 


dN 


n 


dt 


=  NeNn 


i 


a(m\n )  +  z,  P(m|n) 

m>n  m<n 


2> 


+ 


Energy  leve 


Boltzmann 


•  Conserves  number  density 

•  Preserves  energy  in  groups 
through  group 
temperature  description 


Conserved  Variable: 
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Effective  rate: 
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Complexity  Reduction  for  Argon 


Full  Lines  (LANL) 


•  Based  on  quantum  calculations  with 
corrections  for  low  temperature 


NIST  Cutoff 


•  Starts  with  LANL  and  assumes  higher 
excited  states  are  ionized 

•  Cutoff  experimentally  determined 

•  2-3x  reduction 


Electron  Configuration 


•  Groups  based  on  electron  configuration 

•  Uses  uniform  grouping 

•  10-15x  reduction  over  NIST 


Grouping 


•  Boltzmann  or  Uniform  grouping 

•  Saves  20-30%  over  Electron  Splitting 

•  Case  by  case  basis 


Isolated  MUGroups  -B-Electron  Group  - Electron  Group  - Grouping 
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Number  of  Levels  (Runnig  Total) 


A 


Simulation  Setup 


Argon  Level  Grouping  Isothermal  Test 

Setup 


Groupings 


j 


•  Pressure:  4.22  Torr  (5.55xl0'3  atm) 

•  Ion  Temperature:  0.035  eV 

•  Atomic  Density:  1020  1/m3 

•  Ionization  fraction:  10  13 

•  Electron  Temperature:  10  &  100  eV 

•  t  =  [0,106]  seconds 


•  NIST  cutoff  with  electron  grouping 

•  NIST  cutoff  with  electron  grouping 
and  Boltzmann  grouping 

•  NIST  cutoff  with  electron  grouping 
and  Uniform  grouping 
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ln(N/g) 


Argon  Level  Grouping  Isothermal  Test 
Electron  Temperature  10  eV 


Argon  +1 

12  electron  configurations 
9  Isolates  states 
1  group  with  3  states 


Uniform  group  propagates 
errors  to  isolated  states 


24.0  24.5  25.0  25.5  26.0  26.5  27.0  27.5  28.0 

Energy  [eV] 
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ln(N/g) 


Argon  Level  Grouping  Isothermal  Test 
Electron  Jemperatjjre  100  eV 


Argon  +3 

17  electron  configurations 
14  Isolates  states 
1  group  with  3  states 


Uniform  groups  propagates 
errors  to  isolated  states 
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Relative  Difference  (L2) 


Argon  Level  Grouping  Isothermal  Test 

Error 


ODE  Solver  tolerance 
Relative  error  1~4 
Absolute  error  1 
~1200  Radau5  time  steps 
dt  increases  exponentially 


Solution  improves  if 
either  tolerance  is 
decreased  but  at  the 
expense  of  computational 
time.  E.g.  relative  error 
10~6  ->  computational 
time  triples 


Time  [secs] 


Distribution  A  -  Approved  for  public  release;  distribution  is  unlimited;  PA:  #16531 


14 


SL 


Non-Maxwellian  Electrons  CR  -  Ions 

(Preliminary) 


Current  Assumptions 

—  Electron  dominated 
collisions 

—  Single  ionization  level 

—  Isotropic  EEDF 

—  6  function  for  ion 
distribution 

Current  Model  Includes 

—  Elastic  electron  collisions 

—  Excitation/de-excitation 

—  lonization/recombination 


dNn(t) 

dt 
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SL 


Non-Maxwellian  Electrons  CR  -  Electrons 

(Preliminary) 


Current  Assumptions 

—  Electron  dominated 
collisions 

—  Single  ionization  level 

—  Isotropic  EEDF 

—  6  function  for  ion 
distribution 

Current  Model  Includes 

—  Elastic  electron  collisions 

—  Excitation/de-excitation 

—  lonization/recombination 
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Non-Maxwellian  CR  Numerical  Challenges 


•  Long  complex  CR  formulas 


Stiff  equations 
Range  of  scales 
Boundary  conditions 
Multi-dimensional  integrations 


Hydrogen  recombination  example: 

eo  +  Hn  e1  +  e2  +  //+ 

Energy  equation  in  terms  of  the  electron's 
kinetic  energies,  e,  and  ionization  energy  In. 

£o  —  £i  +  e2  +  In 

Evaluating  the  effect  of  recombination  on  a 
single  species  of  hydrogen  requires  the 
evaluation  of  a  2D  integral. 


£o 
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Hybridization  Techniques 


•  Solve  ODE/PDE  with  a  coarse  "C" 
solver  and  fine  "F"  solver 


Simpler  hybrid  (H)  method 


Coarse  propagates  less  information 
and  is  more  computationally 
efficient 

Fine  propagates  more  information 
and  is  more  accurate 


_  ( C[t+dt,t](<?0f needed))  if  h<€ 
[t+dt,t]  |f*[t+dt>t](tf(if  needed))  if  h  >  e 

More  complex  time  parallel  (TP)  method 


User  defined  coarse  error  function 
h(u),  0  <  h  <  1 

when  h=0  coarse  is  accurate 

Q  compression  operator  H2ydti0](um 

R  reconstruction  operator 

h,  £,  Q,  &  R,  are  problem  dependent 


C[/dt,o]Q(“(0))  if  h<e 

(  C[jdt,o] Q(ti{0))  \ 

1  1 


RC[jdt,(j+l)dt\  QF(j +l)dt,jdt]  RC[jdt,Q]  Q  (^(0)) 

;=0 

y=y-i 


RC[jdt,(j+l)dt](-‘(j+i)dt,jdt]  C 

;=0 


/ 


F\jdt,o\(u(0))  if  h>  y 


if  e  <  h  <  y 
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Error  (Scaled) 


Hybridization  Charged  Particle  in  Magnetic  Mirror 


*Submitted  to  Journal 


Computation  Time  (sec) 


Bx  = 


Bz  =  C 


2  Cxz3 


2Cyz3 


a4 


Gyrokinetic  +  Phase  and 
Time  Parallel  method  are 
the  most  efficient 


Simple  blending  could  be 
more  accurate  as  the 
problem  approaches 
stead  state 
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Summary 


9  Collisional  Radiative 

—  Boltzmann  grouping  improves  group  representation  over  applied  uniform  distribution 

—  Minimization  of  error,  which  is  expected  to  accumulate  quickly  during  highly-transient 
durations,  and  acceleration  of  method  makes  Boltzmann  reduction  a  strong  case  for  future 
coupled  simulations 

—  Adaptive  integration  allowed  for  faster  simulations  at  larger  electron  temperatures 

—  Sensitive  to  selected  groups 

—  Grouping  reduces  stiffness 

—  Robust  solver  capable  of  handling,  Te  from  1  -  1,000  eV  and  1016-1024  particles 

—  Initial  work  on  non-Maxwellian  CR  has  begun  and  early  numerical  issues  have  been  addressed 

9  Hybridization 

—  Shows  that  hybridization  technique  can  be  more  computational  efficient  than  the  schemes 
that  comprise  it 
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Future  Work 


*  Collisional  Radiative  Simulations 

—  Further  comparisons  between  reduced  mechanisms  with  QSS 

—  Line  identification  and  width  assignments  in  conjunction  with  experimental  spectra 

—  Apply  Boltzmann  grouping  to  new  CR  Argon  rates  and  test  with  ID  MHD  Argon 
shock  by  Kapper,  et.al.  Future  plans  to  extend  to  2  and  3D  MF 

—  Rerun  previous  LPI  test  case  with  level  grouping;  laser  source  term,  heavy-electron 
elastic  collisions,  and  multi-electron  energy  correction  and  heavy  energy  equation 
(Te,  Th) 

*  Group  selection  through  numerical  analysis,  optimization 

*  Hybridize  Maxwellian  CR  with  QSS 

*  Non-Maxwellian  finite  rate  EEDF  implementation 
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Questions? 
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